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Abstract 

In the present paper we study a singular perturbation problem for a Navier-Stokes-Korteweg model with 
Coriolis force. Namely, we perform the incompressible and fast rotation asymptotics simultaneously, while 
we keep the capillarity coefficient constant in order to capture surface tension effects in the limit. 

We consider here the case of variable rotation axis: we prove the convergence to a linear parabolic-type 
equation with variable coefficients. The proof of the result relies on compensated compactness arguments. 
Besides, we look for minimal regularity assumptions on the variations of the axis. 
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1 Introduction 

In the present paper we continue the study, started in |13| . of singular perturbation problems for 
viscous capillary fluids under the action of fast rotation of the Earth. 

Denoting by p > 0 the density of the fluid and by tt G its velocity field, the mathematical 
model is given by the Navier-Stokes-Korteweg system with Coriolis force 

St dtp -b div {pu) = 0 

St dt{pu) + div ipu (8) n) -b 2^Vn(p) - :;^div ipDu) - z^pVAp -b = 0 , 

Fr Re We Ko 

where the positive real numbers zz and k are respectively the viscosity and capillarity coefficients, 
Du is the viscous stress tensor and pS/Ap is the surface tension tensor; finally, 11 represents the 
pressure of the fluid, which is supposed to be a smooth function of the density only. 
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In system (NSK), the positive parameters St, Fr, Re, We and Ro denote respectively the 
Strouhal, Fronde, Reynolds, Weber and Rossby numbers (see e.g. [8], |16|i: we identify the 
Fronde and the Mach numbers. The term (t{p,u) = ce^ x pu represents the Coriolis operator, 
where c is a suitably smooth function (we refer to Section [2] below for the precise assumptions). 

There are two important features to point out in the mathematical theory of equations (NSK). 
First of all, the capillarity term in the momentum equation gives additional bounds for higher order 
derivatives of the density. Such a property shows up not just in the classical energy inequality, but 
also through the so called BD entropy conservation, a second energy inequality first discovered 
in [1] by Bresch and Desjardins (see also E) for our system, and then generalized by the same 
authors to different models: see e.g. 0 , 0 . It turns out that the BD entropy structure is 
a fundamental ingredient in the theory of compressible fluids with density-dependent viscosity 
coefficients: for instance, we quote here works uni, m, m, m, m- 

On the other hand, one has to remark that the viscosity term is degenerate in regions of 
vacuum, where one loses then any information on the velocity field and its gradient. For this 
reason, for system (NSK) supplemented with a classical barotropic pressure law (this hypothesis 
was however a bit relaxed), in [7] Bresch, Desjardins and Lin proved the global in time existence 
of a modified notion of weak solutions. Namely, under these assumptions stability can be obtained 
only in non-void regions of the physical space, so that one has to test the equations only against 
functions supported in zones where the density does not vanish. This was achieved by (formally) 
using pij: as a new test function, with ?/) £ T>t^x and p the density itself, and passing to the limit 
in the new “weak” formulation. 

However, it turns out that stability can be recovered in presence of additional terms in the 
momentum equation: for instance, some drag forces (like in [3] by Bresch and Desjardins for a 2-D 
shallow water model), or a “cold component” in the pressure law (like in 0 by the same authors, 
where variations of the temperature are taken into account as well); alternatively, it is possible to 
impose some additional integrability conditions on the initial velocity, as done in |26| by Mellet 
and Vasseur. We refer to paper 0 for a complete and interesting discussion on this subjet, as 
well as for some recent develpments (see also [30]). 

Differently to what done in |13| . where the weak formulation of |7| was used, we suppose here 
that the pressure term H is given by the sum of two components, a classical one P{p) = p^ and a 
singular one Pc{p) = for two suitable parameters 7 > 0, 7 c > 0. Therefore (thanks to the 

results of 0 , under much more general assumptions than the ones considered here) one recovers 
existence of global in time weak solutions in the classical sense; let us point out that, alternatively, 
we might have added some drag terms to our equations (as done in 0 ), without substantially 
changing the subsequent analysis. The term Pc is often referred to as cold pressure (see e.g. 
0 ), because it is associated with the zero Kelvin isothermal curve for heat conducting fluids; 
also, singular pressures naturally appears when e.g. Van der Waals type laws are considered. In 
fact, for densities and temperatures close to 0 the properties of the medium drastically change, 
damaging the validity of the equations of motion: the presence of Pc may be seen as a way of 
preserving stability of the model. At the mathematical level, this term gives a control for negative 
powers of the density, which can be used to deduce integrability properties on the velocity field. 
We refer to Subsection 3.1 of 0 for more details and some physical insights. Let us also recall 
other works involving singular pressure laws: for instance, |28| for mixture of fluids with chemical 
reactions, |22| for some lubrication models in one space dimension and the above mentioned work 
0 for compressible Navier-Stokes equations. 

In the sequel, we also assume St = Re = 1 in system (NSK), and we set k = 1 for convenience. 
Moreover, for e £ ]0,1] we take Fr = Ro = e and We = where a £ [0,1]: we end up with 

the system 


{ dtp + div {pu) = 0 

, , 1 1 . . 1 

dt {pu) + div {pu ®u) + Vn(/?) -I- - ^{p, u) - vdW {pDu) - ^ • 
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We are interested in studying the asymptotic behavior of a family of weak solutions (^p^,Us) for 
e —>■ 0. This means that we are performing the incompressible and fast rotation limits simulta¬ 
neously, combining them with effects coming either from vanishing capillarity (for 0 < a < 1) or 
constant capillarity (i.e. a = 0). 

The mathematical study of fluids in fast rotation has now a quite long history, which goes 
back to the pioneering paper [T] by Babin, Mahalov and Nicolaenko for the incompressible Navier- 
Stokes equations. We refer to book |12| and the references therein for a complete analysis of the 
problem for incompressible viscous fluids; in the matter of this, we also quote here papaer [6] and 
|18| . Let us just point out an important aspect in the theory of rotating fluids (see e.g. |12)i: 
the strong Coriolis force has a stabilazing effect on the motion. Namely, in the limit e —0, 
the dynamics becomes constant in the direction parallel to the rotation axis: the fluid moves 
along vertical coloumns (the so called “Taylor-Proudman coloumns”) and the flow is purely 2- 
dimensional, evolving on a plane orthogonal to the rotation axis. We refer also to book |29| for 
some physical background about this problem. 

In the compressible case, various models have been considered, for which the choice of the 
scaling Fr = Ro = e is usually assumed. The classical barotropic Navier-Stokes equations with 
constant viscosity coefficient (not depending on the density, then) were studied in |15) by Feireisl, 
Gallagher and Novotny, and in m by the same authors in collaboration with Gerard-Varet. For 
general ill-prepared initial data, in the former article it was proved the convergence of the model to 
a 2-D quasi-geostrophic equation, by resorting to the spectral analysis of the (constant coefficient) 
singular perturbation operator. In the latter work, instead, still for ill-prepared data, the limit 
system was identified as a linear equation, due to the presence of the centrifugal force having the 
same scaling as the other quantities. Also, in m the singular perturbation operator is no more 
constant coefficients, and the authors had to resort to compensated compactness arguments in 
order to pass to the limit in the non-linear terms. This method was firstly introduced by P.-L. 
Lions and Masmoudi (see |23|-|24|1 in dealing with incompressible limit problems, and borrowed 
in |18| by Gallagher and Saint-Raymond in the context of rotating fluids. 

Goncerning systems having a similar structure as in the above mentioned paper [3] Bresch 
and Desjardins studied the problem for a 2-D viscous shallow-water model with density dependent 
viscosities (we refer also to |19| and m for a similar system, but with constant viscosity coeffi¬ 
cient) and with additional laminar and turbolent friction terms; also, they assumed a vanishing 
capillarity regime. In |21| . instead, Jiingel, Lin and Wu considered more general viscosity and 
capillarity tensors, working in a strong solutions framework, still in space dimension 2 and in the 
vanishing capillarity regime. In both works [3] and [2T], for well-prepared initial data, the authors 
recovered convergence to a quasi-geostrophic equation, by use of the modulated energy method. 

Goming back to the somehow simpler form ([1]), in m we studied the problem in the 3-D 
domain D = M^x ]0,1[ and for general ill-prepared initial data: the improvement was due to the 
use of spectral analysis tools (namely RAGE Theorem and microlocal symmetrization arguments), 
as done in m- We payed attention both to the vanishing and constant capillarity cases: in the 
former instance we recovered the asymptotic result of [3] and |21j . while in the latter we found the 
convergence to a slightly modified 2-D quasi-geostrophic equation: since surface tension effects 
do not vanish in the limit, then they come into play also in the final relation. 

In the present paper we want to continue the previous study, mainly focusing on the case 
of effectively variable rotation axis depending just on the “horizontal variables”, in analogy to 
what done in |18| by Gallagher and Saint-Raymond for the homogeneous incompressible Navier- 
Stokes system. Indeed, when c = 1 (so that C(/9, u) = x pu) the analysis of |13) still applies: the 
presence of the cold pressure is important just for stability of weak solutions, but it does not affect 
the singular pertubation problem (actually, it simplifies things, since it supplies informations for 
u and its gradient). We will come back in a while to the hypotheses for the function c. For 
the moment, let us point out that we restrict our attention to the constant capillarity regime, 
corresponding to a = 0: we are interested here in capturing surface tension effects in the limit; 
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on the other hand, different choices of a can be treated in a similar way. 

For non-constant rotation axis, the singular perturbation operator becomes variable coeffi¬ 
cients, so that spectral analysis tools are no more available. Resorting to the techniques of |18) . 
also used in m for non-constant density prohles, the idea is hence to apply compensated com¬ 
pactness arguments to prove the convergence in the non-linear terms: more precisely, after a 
regularization procedure and integration by parts, we take advantage of the structure of our sys¬ 
tem to find special algebraic cancellations and relations, which enable us to pass to the limit. The 
main novelty here is the presence of an additional non-linear term, due to capillarity: the BD 
entropy structure gives compactness in space for the density, but we miss uniform informations 
in time, so that we cannot pass directly to the limit in it. Nonetheless, it turns out that this item 
exactly cancels out with another one, coming from the analysis of the convective term. Notice 
also that the regularization process itself presents some complications with respect to |18) and 
because we have here less available controls for the velocity field and its gradient. In the end, 
we can prove the convergence to a variable coefficients linear equation (of parabolic type) for the 
limit density, which can be seen as a sort of stream function for the limit velocity field. 

As it was already the case in |18j . the previous arguments work under high regularity assump¬ 
tions on the variation c of the axis. First of all, a control on the gradient of c is needed in the 
computations, in order to use the vorticity equation and to decompose the horizontal part of the 
(approximated) velocity field into the basis given by Vc and its orthogonal V-’-c (recall that c just 
depends on the horizontal variables). Hence, having c € W^’°° seems to be a necessary hypothesis, 
at least for this strategy to work; besides, we will also need to assume that Vc has non-degenerate 
critical points, in a precise sense (as already done in |18|1. 

But this is not all: the regularization procedure creates some remainder terms (essentially, 
commutators between a smoothing operator and the variable coefficient), which we need to be 
small together with their gradient in order to be able to exploit the vorticity formulation. This 
requirement asks for additional regularity on Vc: in the present paper we look for minimal 
smoothness assumptions on it in order to recover the convergence. More precisely, we show 
that it is sufficient for Vc to be //-continuous, for some admissible modulus of continuity /t. The 
proof of this result relies just on fine commutator estimates, which can be obtained, nonetheless, 
in a quite classical way. 

To conclude, we give an overview of the paper. In the next section we set up the problem, 
formulating our working hypotheses, and we state the main results. Section [3] is devoted to 
establish suitable a priori estimates, in the general case of system [NSK)\ besides, this analysis 
allows to justify a technical requirement in |13| for the vanishing capillarity regime. In Section 
|4] we prove our main result, about the singular limit problem in the case of variable rotation 
axis, when assumptions are made on the first variation of Vc. Finally, in the Appendix we recall 
some basic notions from Littlewood-Paley theory (Appendix lA.ip and about admissible moduli 
of continuity (Appendix |A]2]) ; we postpone the proof of some technical results for the BD entropy 
in Appendix I A. 31 

Notations. Let us introduce some notations here. 

We will decompose x £ D := M^x ]0,1[ into x = (x^, x^), with x^ € denoting its horizontal 
component. Analogously, for a vector-field v = (x^, x^) £ we set x^ = (x^, x^), and we define 

the differential operators V/i and div/i as the usual operators, but acting just with respect to x^. 
Finally, we define the operator V^ := {—62 , i9i). 

Moreover, since we will reconduct ourselves to a periodic problem in the x^-variable (see 
Remark 12.II below), we also introduce the following decomposition: for a vector-field X, we write 

(2) Ar(x) = {X){x^) -|- V(x), where {X){x^) := f X{x^,x^)dx^. 

Jt 
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Notice that X has zero vertical average, and therefore we can write X[x) = d^Z{x), with Z 
having zero vertical average as well. We also set Z = T{X) = d^^X. 
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2 Basic definitions and results 

We describe here our main hypotheses and results. 

2.1 General setting 

Let us consider the rescaled Navier-Stokes-Korteweg system 
dtp + div {pu) = 0 

(3) < I II 

dt{pu) + div (pu (g) u) + :^Vn(p) - vdXipDu) - z—pXAp + — (t(p, tt) = 0 
, Fr We Ko 

in M_|_ X O, where O is the infinite slab 


0 := X ]0,1[. 

In the previous system, the scalar function p > 0 represents the density of the fluid, while rt £ 
its velocity field. The quantity i', which we will assume always to be a fixed positive constant, is 
the viscosity coefficient and the operator D denotes the viscous stress tensor, defined by 

Du := ^ (Vm -|- *Vri) . 

Moreover, we denoted by C the Goriolis operator, which takes into account the Earth rotation: 
here, we suppose that <t is given by 

(4) ^{p,u) := c{x^)e^ x pu, 

where = (0, 0 ,1) is the unit vector directed along the x^-coordinate and c is a smooth scalar 
function of the horizontal variables only. Finally, the function n(p) represents the pressure of the 
fluid: for the reasons we explained in the introduction, we suppose here H = P + Pc, where P 
is the classical pressure, given by the Boyle law 

(5) P{p) ■= — P'^ 1 for some 1 < 7 < 2, 

27 

and the second term is the cold pressure component, for which we take the power law 

( 6 ) Pc{p) p”"^" , with 1 < 7 c < 2 . 
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The presence of the 1/2 is just a normalization in order to have 11'(1) = 1: this fact simplifies 
some computations in the sequel. 

In view of the analysis of the singular perturbation problem, we supplement system ([3]) by 
complete slip boundary conditions, in order to avoid the appearing of boundary layers effects. 
Namely, if we denote by n the unitary outward normal to the boundary dO, of the domain (simply, 
on = = 0} U {x^ = 1}), we impose 

(7) iu-n)\Q^ = = 0, {Vp-n)\Q^ = dspian = 0, {{Du)n x n) = 0. 

Remark 2.1. Equations Q, supplemented by complete slip boundary boundary conditions, can 
be recasted as a periodic problem with respect to the vertical variable, in the new domain 

X with := [-l,l]/~, 

where ~ denotes the equivalence relation which identifies —1 and 1. Indeed, the equations are 
invariant if we extend p and as even functions with respect to x^, and as an odd function. 

In what follows, we will always assume that such modifications have been performed on the 
initial data, and that the respective solutions keep the same symmetry properties. 


2.2 Existence of weak solutions 


Here we will always consider initial data {pq^uq) such that po ^ 0 and 


( 8 ) 


l(po - 1) e L^{n) 


, V y/pQ 


1 / 1 


VWe 


and —(—-lie 
Fr Vpo 

Vpo e l2(0). 


At this point, let us also introduce the internal energy functions h{p) and hc{p), such that 

P'(n) 

h''{p) = - = p'^~‘^ and h{l) = h'[1) = 0, 

P'(o) 

h"{p) = — - = p~'^‘^~‘^ and ^c(l) = ^c(l) = Oj 

P 

and let us define the classical energy 

(9) E[p,u\{t) := j^(^-^h{p) + -^hc{p) + \ 

and the BD entropy function 


( 10 ) 


F[p]{t) 




Y^y/P? dx. 


Finally, let us denote by E[pq,uq] = E[p,u]{0) and by E[pq] = T’[p](0) the same energies, when 
computed on the initial data (poi^^o)- 

We now give the definition of weak solution to our system. The integrability properties we 
require (see points (i) and (ii) below), as well as conditions ([8]) for the initial data, will be justified 
by energy estimates, which we will establish in Subsection 13.11 


Definition 2.2. Fix initial data {po,uo) which satisfy the conditions in ([8]), with po ^ 0. 

We say that (p, u) is a weak solution to system ([3])-([7]) in [0,T[xn (for some T > 0) with 
initial datum (po,Mo) if the following conditions are verified: 
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(i) /9 > 0 almost everywhere, and one has the properties Fr — l) E ([0, T[; L'^(n)), 

Fr-i(l/p-l) G L°°{[0,T[;L^-{n)), We-^/^Vp and G L°°{[0,T[; L^{Q)) and 

We-^/2y2^ g L2([0,r[;L2(O)); 

(ii) G L°°([0,r[;L2(O)) and G ([0, r[; ; 

(iii) the mass and momentum equations are satisfied in the weak sense: for any scalar function 
4> G P([0, r[ xff) one has the equality 


/7 


pdtcj) + pu ■ V(f> ] dxdt = 


/ po</>(0)dx, 
Jn 


and for any vector-field il> G P([0,r[xn;one has 

(11) [ PoUq ■ ip{0) dx = I I ( —pu ■ dtip — pu®u: Vip — —^ n(p) divip + 
Jn Jo Jn\ Fr^ 


Fr^ 

1 1 c(x^') 

+ vpDu: Vijj + pApdivil) -h ApVp ■ ijj H--— 

We We Ro 


X pu ■ 'll)] dx dt. 


For such initial data, the existence of weak solutions to system Q is guaranteed for any fixed 
value of the positive parameters Fr, We and Ro. 

Theorem 2.3. Let 7 c = 2 m ([6|) and c G in (|4]). Fix the value of the Froude, lFe6er 

and Rossby numbers, and consider an initial datum {pQ,uo) satisfying conditions ([S]), with po > 0. 

Then, there exits a global in time weak solution {p,u) to system Q, related to {pq,uq). 

The previous result can be established arguing exactly as in [9] , so we omit its proof. Actually, 
the result of [9] holds true under more general assumptions than ours (as for the cold component 
of the pressure and the viscosity coefficient, for instance). 

Remark 2.4. • The hypothesis 7 c = 2 is assumed just for simplicity here, but, as remarked 

above, it is not really necessary for existence. 

• The condition c G W^’°° is important in order to take advantage of the BD entropy structure 
of our system, see Paragraph 13.1.21 However, it can be deeply relaxed at this level: see also 
the discussion at the beginning of Subsection 14.41 


2.3 The singular perturbation problem 

We get now interested in a singular perturbation problem for system ([3]). Namely, we want to 
study the incompressible and high rotation limit simultaneously, both in the regimes of constant 
and vanishing capillarity (in the same spirit of the analysis of m)- 

For doing this, we consider a small parameter £ G ]0,1]: we set Fr = Ro = e, We = 
for some 0 < a < 1. Notice that the constant capillarity regime corresponds to the choice a = 0, 
while in the other cases we are letting also the capillarity coefficient go to 0. 

Therefore, we end up with the equations 


r dtp + div (pu) = 0 

( 12 ) < 11 1 

ydt (pu) + div [pu (8) Ii) + ^ Vn(/9) -k - C(/9, u) - z^div (pDu) - ^ > 

Here we will consider the general instance of ill-prepared initial data (p, tz)|j_Q = (po,£)^o,£)- 
Namely, we will suppose the following assumptions on the family (po,£ ; '*^o,£)c>o’ 
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(i) /9o,e = 1 + e?’o,£) with (j’o.e)^ C H^{Q,) r\L°°{Q,) bounded; 

(ii) l/po,e = 1 + eao,£, with (ao,e)^ C bounded; 

(hi) (iio,e)j C L'^iyt) bounded. 

Up to extraction of a subsequence, we can suppose that 

(13) ro,£ ^ ro in H^{9) , ao,£ ^ ao = - in U^(ll) , 'Uo,£ ^ uq in L^(r2), 


where we denoted by ^ the weak convergence in the respective spaces. 


Remark 2.5. The property oq = —ro immediately follows from the weak convergence. Indeed, 
for any test function (p S by definition of ao,£ we have 


/ 

Jn 


ao,£ p dx 


1 

e 



1) ipdx 



ro,£ p dx 



ro,e ao,e ^ dx . 


The left-hand side of the previous equality converges to aopdx; as for the right-hand side, the 
former term converges to — ropdx, and the latter goes to 0 for e —0. 

Remark 2.6. Notice that the choice of constant density profile in the limit, i.e. p = 1, is 
consistent with the balance of forces acting on the system. As a matter of fact, p is identified as 
a solution of the static problem 

Vn(p) = -e^^pVAp, 

and this relation implies, up to an additive constant, 

— e^“Ap = H(p), with H(p) := J [H' (a)/a) da . 

Of course, p = 1 solves the previous elliptic equation. However, depending on the non-linearity 
n, when a = 0 one is led to consider also non-constant limit density profiles: this was already the 
case in m for a barotropic Navier-Stokes system with Coriolis term, when the centrifugal force 
is assumed of the same order as the pressure and the rotation. 

The analogue for capillary fluids will be matter of future studies. 


Furthermore, we need to slightly modify the definition of weak solutions: namely, in addition 
to the conditions of Definition 12.21 we also demand that the weak solutions are constructed in 
such a way to satisfy relevant uniform bounds in e. More precisely, we set E^[p,u] and ^^[p] the 
energies defined in (l9|) and (fTOl) respectively, where we take the scaling Fr = Ro = e, We = : 

hence we require that, for almost every t £]0, T], the following inequalities hold true: 


(14) Es[p,u]{t) + [ f P iDuf' dxdr 

Jo Jn 

(15) Ee[p]{t) + ^ J^P'{p)\V^f dxdr -h 


< 


< 


Ee[PO,Uo] 

C{l + T), 


where the constant C depends just on the triplet (^^^[po, uq], ^^[po], i^). 

From now on, we will focus only on the regime of constant capillarity, i.e. a = 0. Indeed, our 
main goal is to capture the effects of the surface tension in the asymptotics, which seems to be a 
new feature in this kind of studies. However, we remark that the vanishing capillarity limit (i.e. 
a g] 0, 1]) can be dealt with as in |13| when c = 1, or by similar arguments as in Section 0] when 
c is non-constant. 

First of all, let us consider the case of constant rotation axis, namely when c = 1 in (j3|). 



Theorem 2.7. Let 1 < 7 < 2 m ([5]), a = 0 and ^(p,u) = x pu in (1121) . 

Let (po^e, he initial data satisfying the hypotheses (i) — (ii) — (iii) and d, and let 

(pe , tie) he a family of corresponding weak solutions to system (I12p -(j7|l in [0,T] x 12, in the 
sense of Definition \2.S[ Suppose that inequalities dH-dED hold true and that the symmetriy 
properties of Remark \2. 1\ are verified. We also define r^ := {p^ — 1). 

Then, up to the extraction of a subsequence, , Ug)^ weakly converges (in suitable en¬ 
ergy to {r,u), where r = r{x^) and u = (ti^(x^), 0 ) are linked by the relation = 

(Id — A/i) r. Moreover, r is a weak solution of the modified quasi-geostrophic equation 

(16) 5i(ld - A/, + A|)r + V^(ld -A;,)r • Vu^lr + ^Ai(ld -Afc)r = 0 

supplemented with the initial condition r|^=o = ^ 0 ; where Fq G is identified by 


(Id - Ah + Al) ro = f (wg + ro) dx^ . 

Jo 

We omit the proof of this result here, because it goes along the lines of the one given in Ha¬ 
lts main ingredients are the spectral analysis of the (constant coefficient) singnlar pertnrbation 
operator and an application of the RAGE Theorem. We remark that the RAGE Theorem allows 
to deduce strong convergence properties in suitable norms, and this is the key to pass to the limit 
in the non-linear terms. 

Let us now consider the case of effectively variable rotation axis. Eor technical reason, anal¬ 
ogously to what done in |18) . we need to assnme that the fnnction c has non-degenerate critical 
points: namely, we will suppose 


(17) 


lim C({x^ G 
5^-0 V I 


|Vfcc(x'^)| < j}) = 0, 


where we denoted by C{0) the 2-dimensional Lebesgue measure of a set O C M^. 

Also, we will require that the gradient of c is /i-continuous, for some admissible modulus of 
continuity p: we will recall the precise definition in Appendix IA.2I 
For notation convenience, let us also introduce the operator 

S),(/) := Z2,(c-iV^/) = ^ (V;. + *V,) (w^ V^/) 
for any scalar function / = f{x^). 

Theorem 2.8. Let 1 < 7 < 2 OT 0 = 0 and <t{p,u) = c{x^)e^ x pu in (fT^ . where 
c G iy^’°°(M^) is 7 ^ 0 almost everywhere and it verifies the non-degeneraey condition m- Let us 
also assume that V/jC G C^(M^), for some admissible modulus of continuity p. 

Let (/ 0 o,£)^ o,£)£ be initial data satisfying the hypotheses (i) — (ii) — (iii) and (fT^ . and let 
{pe, Ue)^ be a family of corresponding weak solutions to system ([12])-([7]) in [0, T] X Ll, in the 
sense of Definition 12.21 Suppose that inequalities dH-dED hold true and that the symmetriy 
properties of R.emark \2. 1\ are verified. Define r^ '.= {pe — 1) os above. 

Then, up to the extraction of a subsequence, one has the following convergence properties: 

(a) re ^ r zn ([ 0 , T]; F^f^)) G LWf),T]-H‘^{D)), 

(b) y/pfue u in L°° [[0,T]-, L‘^{LI)) and^fpfDue Hu m L^([0, T]; L^(f2)), 

^See points (a) and (b) in Theorem 12.81 
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where, this time, r = r(x^) and u = (m^(x^),0) verify the relation c{x^)u^ = (Id — A/j)r. 
Moreover, r solves (in the weak sense) the equation 



supplemented with the initial condition r|^=o = ''’0; where tq is defined by 

fo - V/i(ld - A/i)ro^ = J (^curl/,(c“^ Uq) + dx^ 

Notice that we have the identity 


o D,(/) = VM 


where we used the notations div / and V • / in an equivalent way. We also remark here that, for 
c = 1, this operator reduces to (1/2)A^/. 

Remark 2.9. As already pointed out in |18) (see also |14)h the limit equation is linear in the case 
of variable rotation axis. Indeed, the limit motion is much more constrained in this situation, and 
correspondingly the kernel of the singular perturbation operator is smaller than in the instance 
of constant axis. 

Also, notice that having non-constant c makes variable coefficients appear in the limit equation. 


3 A priori estimates 

The present section is devoted to show uniform bounds for smooth solutions of our system, 
written in the general form ([3]). Throughout this section we will suppose Fr, We and Ro to be 
hxed. However, we will keep track of them in the first paragraph: this is relevant in view of the 
analysis of Section 01 

Besides, the next computations will give complete justification to a technical assumption in 
m (see Remark 13.81 belowl. 

3.1 Energy estimates 

Suppose that (p, u) is a smooth solution to system ([3]) in M_|_ X fl, related to the smooth initial 
datum (po)'^^o)- We establish here energy estimates for {p,u). 

Remark 3.1. The construction given in e.g. [9] ensures the existence of smooth approximated 
solutions, which converge to a weak solution of our original system and which are compatible with 
the uniform bounds given by the BD entropy structure of the equations. 

Therefore, the estimates established here will be inherited by the family of weak solutions we 
are going to consider, see Subsection 14.11 

3.1.1 Classical energy 

First of all, let us show the classical energy conservation. 

Proposition 3.2. Let {p,u) he a smooth solution to system Q in M_|_ x 12, with initial datum 
{pOiUq). Then, for all t S M_|_, one has 

d f 

—E[p,u]{t) + j p\Du\^ dx = 0. 

dt Jo 
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Proof. First of all, we multiply the second relation in system (|3|) by u: by use of the mass equation 
and the fact that (L{p, u) is orthogonal to u, we arrive at the equality 


^ ^ j^{P'ip) + Pcip))'^P-'^dx + ^j^pDu : Vudx = 0. 

On the one hand, we have the identity Du : Vrt = \Du\‘^; on the other hand, multiplying the 
equation for p by h'{p)/¥r‘^ gives 


1 


1 d 


o I P'(p)Vp-udx = —TT — I h(p)dx. 
Fr2 Jn ^ Fr2 dt 


Notice that an analogous equality holds true also for the cold part Pc of the pressure. 

Putting the obtained relations into the previous one concludes the proof of the proposition. □ 


From the previous energy conservation, we infer the following bounds. 

Corollary 3.3. Let 7 c = 2. Let {p,u) be a smooth solution to system Q in M_|_ X Ll, with initial 
datum (pO)^o)j assume that E[pQ,uo] < +oo. Then one has the following properties: 

y/pu, -^Vp e L°°(M+;L2(0)) and ^Du e L^{R+-L^{n)). 

V We 

Moreover, we also get 

- 1) e L°°(]R+;LT'(f2)) and ^ ^ L°°{W+-, L^{M)) . 

Proof. The first properties are immediate, after integrating the relation of Proposition 13.21 in time. 
So, let us focus on the density term. 

For 1 < 7 < 2, Lemma 2 of |21j . combined with Proposition 13.21 tells us 


IIP-llH-IL-.) <C{Pr’ + Fr^), 

which immediately gives the first property. 

Let us consider negative powers of the density: for 7 c = 2, it is easy to see that 


Hc{p) ■■= 7c (Tc + 1) hc{p) 


1 

P 



> 0 


for all p > 0, and this concludes the proof of the corollary. 


□ 


3.1.2 BD entropy 

We want now to take advantage of the BD entropy structure of our system. Let us start with a 
lemma, which is the analogue of Proposition 3.3 in m- For the sake of completeness, we give the 
complete proof in Appendix lA.31 


Lemma 3.4. Let {p,u) be a smooth solution to system ([3]) in M_|_ X Ll, related to the initial datum 
{pOiUoj. Then there exists a “universal” constant C > 0 such that, for all t € 1R+; one has 



p |tt + i^VlogpF dx 


+ 




2 dxdr + ^ 


0 Jn 


n^(p) IV-^pI^ dxdr < 


< C[F[pq]-L E[pQ,U[f\) + — 


c(x^) X u -Vpdx dr 


0 Jn 
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Let us point out here that the last term on the left-hand side of the previous relation can be 
also written as 

(20) ^ l^{P\p)+P^ip))\V^\^ = ^^(p7-2 + ^->-2)|Vp|2 


C.y V 




for some positive constants C-y and In particular, in our case 7 c = 2, = 1. 

Next, let us give now some estimates for the density (again, see the proof in Appendix IA.3p . 

Lemma 3.5. There exists a “universal” constant C > 0 such that 

\\p - 1|Il-(l 2) < C (^Fr + (1 - 12(7)) , 

where we have set 12 ( 7 ) = 1 if 7 = (md 12 ( 7 ) = 0 otherwise. 

Moreover, for any 0 < 5 < 1/2 and any 1 < p < 4/(1 + 26) one has 

11'’“- Cp((Er+(1 - ll 2 ( 7 ))VWe) t + 

where we have defined q := 4/((l + 26)p') £ [1,4/(1 + 25)], and the constant Cp depends also on 
the value of p. 

We are now ready to control the Coriolis term in the estimates of Lemma 13.41 We will give 
two different types of estimates: the first one is useful in order to justify Definition 12.21 of weak 
solutions. The second kind of inequalities, instead, will be relevant in studying the singular limit 
problem (see Section |3|). Again, the proof is postponed to Appendix I A. 31 

Lemma 3.6. (i) There exists a positive constant C, just depending on E[pq,uq], such that 


ix 

Ro 


/'/ 


x^) e^ X u -Vpdx dr 


< 


C_ 

Ro 


Jo 


1/2 


dr. 


(a) Moreover, one has the following estimate: for any 1 < 7 < 2, 


v 

Ro 


/7 c(x") 
Jo Jn 


e^ X u ■ Vpdxdr 


< c.(i + i)h±^('i + FV + x/my'+ 


Ro 


+ C'z/(l + t)(We)3/5 ( ^ 3^ 


Alternatively, in the particular instance 7 = 2 , one can get the different bound 


V 

Ro 


/Fr 

8/5 




n c(x^) e^ X u - Vpdx dr 
i 

+ C'z.(l + t)(We)3/5 W^PWhfiL^) + P^pfLl(L^) ■ 


By the previous lemmas, we easily infer the BD entropy estimate for our system. 

Proposition 3.7. Let {p,u) be a smooth solution to system ([3]) in x D, related to the initial 
datum (po; • 
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(i) There exists a constant C > 0, just depending on E[po,uo] and -F[/9o]; such that, for all 
T £ M_|_ fixed, one has 


supF[/ 9 ] + 

[0,T] 


We 


fl 


2 , ,, ^ 
dxdt + —^ 

JV2 



0 Jn 


n'(p) dxdt < C + C 


T Y 


(a) Alternatively, for any 1 < 7 < 2 and any T G fixed, we have the estimate 

supF[p] + ^ [ [ |VV|^ dxdt + [ [ n'(p) \Vy/p\^ dxdt < 

[o,T] We Jo Jn Fr Jo Jn 

< C + Ciu{l + T)Q + C2 v{1 + T) 0^/5 , 

where 0 = 0(Fr,Ro,We) = ^Fr + VWe^ / Ro and the positive constants Ci and C 2 are 
uniformly bounded for (Fr,We) varying in some compact set [0,Fr] x [0,We]. 

In the case 7 = 2 , one can derive also an analogous inequality, with the right-hand side 
replaced by the quantity 

C + CvTQ^ + CivVtQ + C 2 l/(l + T) 0 ®/^ 

where Ci and C 2 have the same property as above, and we have defined 0 := Fr/Ro. 


Proof. Our starting point is inequality (|19p : we estimate from below the first term in the left-hand 
side by 

^ [ p{t) + z^Vlog/9(t)p dx > F[p]{t) -]- ! p{t)\u{t)\^ dx, 

^ Jn ^ Jn 

where the last term can be clearly moved on the right-hand side of fll9p and controlled by E[pq,uq\ 
(recall Proposition I3.2F 

Let us now consider the rotating term, and bound it by the first inequality in Lemma I3.H1 
Notice that, at this point, one could use Young inequality to get < 1 -|- R and then apply 
Gronwall lemma: this would give an exponential growth in time for E and the other terms we 
want to control. So we prefer to argue in a hner way. 

More precisely, for any t G [0, T] we have the bound 


_C 

Ro 



Ct 

< — 
Ro 


(^supF[p]^ < 

\m J 


Ct 


Ro 


2 + ^ supF[p]. 

^ [0,i] 


and then, from (1191) . we hnd 


Taking first the supjo^r] of fFe member on the right, and then the supp of the member on the 
left, we deduce item (i) of our statement. 

On the other hand, if we bound the Coriolis term in (I19p using the second type of estimates 
of Lemma [3.61 for any 1 < 7 < 2 we easily find the former inequality in item (ii). 

In the special case 7 = 2, we can alternatively use the last bound of Lemma 13.61 since 
P'(p)|Vy^P = |V/ 9 p, we can absorbe the last two terms of the right-hand side into the left-hand 
term. Then, the inequality with 0 immediately follows. □ 


Remark 3.8. The previous proposition suggests to take 7 = 2 when considering the vanishing 
capillarity regime. As a matter of fact, if We = for some 0 < a < 1, uniform bounds in e 

seem to be out of reach without resorting to 0 rather than to 0 . 
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From Proposition 13-Ti we deduce the next statement. Notice that, here below, the last assertion 
derives also from Lemma [3.51 

Corollary 3.9. Let (/?, u) be a smooth solution to system Q in M_|_ x Ll, with initial datum 
{po^uo'j, and assume that E[pQ,uo] and F’[/Oo] are finite. Then one has the following bounds: 


'vvp e LZ{^+-L\n)) 



In particular, the quantity (p — l) /VWe belongs to . 

3.2 Additional bounds 

In the present paragraph we show further properties, which can be deduced from the previous 
controls of Corollary 13.31 and Corollary 13.91 

In the sequel, we will still keep track of the dependence of the various quantities on the Froude, 
Rossby and Weber numbers. On the contrary, we will drop out the dependence on u, since for us 
it will be always a fixed positive parameter. 

Hence, let us fix a T G M_|_. First of all, from Corollaries 13.31 and 13.91 above, we get 
IQ - l) G Lf{L\Ll)) n Ll{H\a)). 

Then, since \ l/yfp — l|<|l//? — l|,by Sobolev embeddings we also infer 

- i) £ L^(L\a)) r, L^tiEa)) ■ 

Therefore, thanks to energy estimates of Proposition 13.21 we deduce first that 
(21) u = y/pu + '/P'^ £ ’ 

and hence, by an analogous decomposition, also that 

y/pDu G n (^Ll{L^ + L^)y 

Notice that, in particular, Du belongs to therefore, Sobolev embeddings implies the 

additional property u G 

We turn now our attention to the quantity pu. Exploiting the same decomposition as above, 
together with the control on p — 1 provided by the BD entropy conservation (recall 

Corollary 13.9p . we have 

(23) pu = y/pu + {y/p — 1) y/pu G (L^ + C (L^) +T^(L^)^ . 

In particular, this implies that pu belongs also to C L|n(L^). On the other hand, for 

its gradient we have 

D{pu) = pDu + uDp = yfpDu + {y/p- 1) yfpDu + yfpuDyfp. 

The first two terms are in (L^ + L^/^) ^ while the last one is in (L^) . Then we find D{pu) G 

(4{l= + l=/^) + l??(l‘)) 


(22) Du = y/pDu + 


Vp 


- 1 
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Finally, let us consider the quantity u: arguing exactly as in (1231) . we get 

(24) + (p - 1) G L??(l2 + l 3/2) n +L^(l2)) . 

Furthermore, we notice that 

(25) D{p^^^u) = p^Du + ^^/puDp 

= -s/J) Du + (p — 1) -s/J) Du + — ■\/pu Dp G Ly (ii^(n) + (^)) ■ 

Indeed, the first term in the right-hand side of the last relation clearly belongs to L|,(L^); more¬ 
over, by Corollaries 13.31 and 13.91 and Sobolev embeddings, the second and the third terms are in 
L|, (L^/^) . Therefore, by use of Proposition I A.4[ we deduce that u belongs also to L|, (L^(n)). 

4 The singular perturbation problem 

In the present section we prove Theorem 12.81 We first establish uniform bounds for the family 
(pe, Ug) of weak solutions we are considering. Then, we study the singular pertubation operator, 
showing constraints on the limit-points of this family. Finally, we pass to the limit in the weak 
formulation of the equations, proving the convergence to equation (|18p . 

4.1 Uniform estimates 

By use of the analysis of Section [3] (keep in mind also what said in Remark 13.ip , we establish 
here uniform bounds for the family of weak solutions (p^jUg)^. Recall that we have fixed 7 c = 2, 
0 = 0 and 1 < 7 < 2 . 

First of all, we consider the energies E^\p,u\ and F^\p], defined respectively by (II4p and (|15p . 
We remark that, under our hypotheses on the initial data, we deduce the existence of a “universal 
constant” Cq > 0 such that 

Ee[Pe,Ue]{^) + Fe[pe]iO) < Cq . 

Then, by use of Corollary 13.31 we immediately infer the following properties. 

Proposition 4.1. Let (pcUg)^ be the family of weak solutions to system m considered in 
Theorem \2.8l Then it satisfies the following bounds, uniformly in e: 

G L°°{R+-,L\n)) and Du^ G L^R+-L^Q)) 

for the velocity fields, and for the densities 

-(pc - 1) G L-(R+;L^(0)) , i (1 - l") E {R+; L\n)) , 
ivpc G L“(M+;L2(5d)) . 

Remark 4.2. In particular, under our assumptions we always have 

\\Pe - l|ll,°o(M+;L2(0)) < Ce. 

As for the BD entropy structure. Corollary 13.91 implies the following estimates; as for the last 
sentence, one has to use Lemma lA.3l in the appendix. 
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Proposition 4.3. Let be the family of weak solutions to system m considered in 

Theorem \2.8\ . Then one has the following bounds, uniformly for £ > 0; 

'ivV,. \^{pT) , 6 LL(E+;i^(Si)). 

In particular, the family (e“^ {p^ — 1))^ is bounded in for any 2 < p < 4. 

Moreover, arguing exactly as in Subsection 13.21 we can establish also the following bounds, 
uniformly in e. First of all, by the decompositions (|2ip and (1221) we immediately have 

{ue)^ C 

[Due)^ c (4 (l2)+4(l 3/2)) n + . 

In particular, is uniformly bounded in therefore, by Sobolev embeddings we 

gather also the additional continuous inclusion {ue)^ C 
Furthermore, we also infer the uniform bounds 

{p,u,)^ c l??(l2 + l 3/2) n (l??(l2)+4(l2)) 

(26) (^D{p^Ue))^ C • 

In particular, we deduce that is a bounded family in L^(L^^^) n 

For the sake of completeness let us also establish uniform bounds on quantities related to 
Ue- First of all, arguing exactly as in (1241) . we get 

{pTuf)^ C L??(L 2 + L 3 / 2 ) n (l??(L 2 )+ 4 (l 2 )); 

on the other hand, analogously to ([2^ . we have also the uniform embedding 

(27) {D{pll^ue))^ C 4(l2(F!) + l3/2(0)). 

Therefore, by Proposition I A. 41 we infer that is uniformly bounded in L|,(L^(n)). 

Finally, from this fact combined with the usual decomposition y/flf = 1 + {y/Pe — 1) and 
Sobolev embeddings, it follows also that 

(p^u,)^ C Ll{L\n)). 

4.2 Study of the singular perturbation operator 

In the present subsection we identify some properties and constraints on the limit points of the 
family of weak solutions (pg,^^)^. 

First of all, by uniform bounds we immediately deduce that p^ —)• 1 (strong convergence) in 
i7^(n)) n i7^(n)), with convergence rate of order e. So, we can write p^ = 

1 + erg, with the family (r^)^ bounded in the previous spaces. Then we infer that 

(28) r, - r in L“(M+; 77^^^)) n lL(M+; ^^(F!)) . 

In the same way, if we dehne := (l/pe — l)/e, we gather that (ue)^ is uniformly bounded 
in L-(M+;L 2 ( 0 )) n LL(M+ ; 77^(12)). So it weakly converges to some a in this space: arguing 
as done in Remark 12.51 it is easy to check that 

(29) ae ^ -r in L°°(M+; l 2(0)) n Lf„,(M+; 77^(12)) . 


16 





Again by uniform bounds, we also deduce 


(30) in lL(M+; 4 /^( 0 )) 

and Dus Du in L^^^(]R+; , where we have identified (L^)* with L°°. 

Notice also that, by uniqueness of the limit, we have the additional properties 


^/pPUe ^ U 

in L°°(M+;L2(L!)) 

PeUe ^ U 

in LL(M+;L2(L!)) 

^/pfDue Du 

in l2(M+;L2(H)) , 


where ^ denotes the weak-* convergence in L°° (M_|_; L^(n)). 

Let us find now some constraints the limit points (r, u) have to satisfy: the following result 
can be seen as the analogue of the Taylor-Proudman theorem in our setting. 

Proposition 4.4. Let (^p^,Us)^ be a family of weak solutions to system (I12p -([71l. each one related 
to the initial datum (/Oo ,£)and fulfilling the hypotheses fixed in Section\^ 

Let us define := £~^ {p^ — 1), and let {r,u) be a limit point of the sequence (rs,Us)^. 

Then r = r(x^) and u = with div^u^ = 0. Moreover, we infer the properties 


cu^ = V^(ld — Affjr and ■ VhC = 0. 

Proof. First of all, we test the mass equation on any (j) S iD([0,T[ xfi): using the decomposition 
Pe = 1 + eJ'E, we get 


— e 




redt(t> 



PeUe • V 0 


£ / ro^e 4>{0) . 

JQ 


From this relation, letting e —)• 0, by uniform bounds and convergence properties established 
above, we deduce that f^u ■ V(f = 0, which in turn implies 


(31) 


div u = 0 almost everywhere in [0,T] x 0. 


Let us consider now the velocity field: for any G P([0, T[ xO; , we use e ip as a test- 
function in the momentum equation (|lll) . By use of uniform bounds, we see that the only integrals 
which do not go to 0 are the ones involving the pressure 11, the Coriolis force Gl and the capillarity 
term: let us focus on them. 

We start by dealing with the classical part P of the pressure: we can write 

-[ [ VP{ps)-f; = -- [ [ (Pip,)-Pil)-P'il)ips-I)]divf; + -P'il) [ [Vps-i,. 
^ Jo Jn ^ Jo Jn\ J £ Jo Jn 

The quantity P{pe) — P(l) — P'{^) {Pe ~ 1) can be controlled by the internal energy h{p^): since 
h{p^)/e‘^ is uniformly bounded in L^(L^) (see Propositions 13.21 and HU), the first integral tends 
to 0 for e ^ 0. Hence, thanks also to (1281) we find 

- r [ vp{pe)-fj p'ii) r [ Vr-V’. 

£ Jo Jn Jo Jn 

An analogous decomposition allows us to treat also the cold part of the pressure: we find that 
(1/e) VPcips) ■ Ip converges to P^il) /q^ Jq Vr • ip. Therefore, since P'(l) = Pc(l) = l/^i in 

the end we infer the convergence, in the limit for e going to 0, 

- r [ vn(p,) -p, ^ r f vr-p,. 

£ do Jn Jo Jn 
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As for the rotation term, uniform bounds immediately imply 



^{Pe,Ue) ■ V' 




ce^ X u ■ 'Ip . 


Finally, let us consider the integrals coming from the capillarity tensor. The former one can 
be treated writing, as usual, = 1 + (p^ — 1 ): therefore we get 


1 

e 



PeApe ■ div'Ip 


fl 


Ar ■ div 'Ip , 


where we used the uniform bound on Ap^ and the fact that — 1 is of order e in e.g. 

On the other hand, since also Vpe is of order e in one easily gets 


- f [ ApeVpe ■ 'Ip —0. 

^ Jo Jn 

Let us sum up all these informations: from the momentum equation, when e ^ 0, we have 
found that the limit point (r, u) has to verify the relation 

(32) c(x^)e^ X n + V(ld — A)r = 0, 

which means in particular 

(9i(ld — A)r = cu^ 

< (92 (id — A)r = —c'u^ 

_ 83 (id — a) r = 0 . 

This relation immediately implies that (id — A)r = (id — A)r(x^) depends just on the horizontal 
variables. Hence, d^r fulfills the elliptic equation —Ad^r + d^r = 0 in H: by passing in Fourier 
variables on x T^, we deduce that 


(33) d^r = 0 r = r(x^). 

On the other hand, by (|32p we get that also cvP' has to depend just on the horizontal variables; 
taking the 83 derivative (since c = c(x^) and c 7 ^ 0 almost everywhere on M^) implies that 

yh _ yh^^hy 

From this property and the divergence-free condition (|3ip . it follows also that 8 ^u^ = 0, i.e. 
83 U^ = 83 U^{x^). On the other hand, by periodicity Jrj,i 83 VP' dx^ = 0, which entails 83 VP = 0 
by decomposition ([2]), and then = vP‘{x^). By use of the complete-slip boundary conditions 
in the end we deduce that vP = 0 , which in turn gives 

(34) u = (u^(x^),0) , with div/^tt^ = 0. 

Finally, let us apply the rot operator to equation ([32]): we obtain 

(35) divh[c{x^) u^) = 0 u^-VhC = 0, 


where us have used the just proved property (1341) . □ 

Remark 4.5. Notice that, by the previous proposition and the fact that r, u £ (]R_|_; L^(H)), 

we get that actually r £ L°°(R+;R^(H)). 
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4.3 Passing to the limit in the weak formulation 

By Proposition 14.41 we can define the singular perturbation operator 
A : L^(n) X L^(fl) X 

(r,V) ^ (^divP , c(x^)e3 X P + V(ld - A)r) , 


which has variable coefficients: so, spectral analysis tools (employed in |13| ) are out of use here. 
Then, in order to prove convergence in the weak formulation of our equations, we will resort then 
to a compensated compactness argument. 

This technique goes back to works by P.-L. Lions and Masmoudi about incompressible limit 
(see e.g. |23)-|24jf : it was introduced for the first time in the context of highly rotating fluids 
by Gallagher and Saint-Raymond in |18) . There, they dealt with incompressible Navier-Stokes 
equations with variable rotation axis. With the same strategy, in paper m Feireisl, Gallagher, 
Gerard-Varet and Novotny studied the case of a compressible Navier-Stokes system with Earth 
rotation and centrifugal force, when the fixed limit density profile is supposed non-constant. 

Most of our analysis will follow the one performed in m-m Notice however that here we 
have to deal with an additional term, coming from the presence of capillarity. 

Let us consider tests functions cj) E P([0,r[xf4) and ij: G P([0, T[ xf2;such that the 
couple {4>,'4>) belongs to Ker.A. Recall that, by Proposition 14.4[ they satisfy 


div 'tp = 0 and c(x^) x V’ + V (Id — A) (/> = 0 . 

In particular, ip = [ip^, O) and (p just depend on the horizontal variable G and they are linked 
by the relation cip^ = V^(ld — AfPjcp. Finally, combining this property with the divergence-free 
condition for ip, we infer also that V^(ld — Ah)(p ■ V^c = 0. 

First of all, we evaluate the momentum equation on such a ip: taking into account the previous 
properties, we end up with 



Po,6 uo,e ■ ^p{0)dx 



(^-Pe Ue ■ dtip 


PeUe®Ue : Vip 


+ 1^ Pe DUe :\/lp + ^ Aps Vpe ■ Ip + 



X PeUe 



dx dt. 


Notice that the dt and viscosity terms do not present any difficulty in passing to the limit. On 
the other hand, the rotation term can be handled by use of the weak form of the mass equation, 
tested on 0 = (id — Ah)(p: we get 


1 

e 



X PeUe - Ip 


~7 [ [ 

£ Jo Jn 

-[ [peu'p-Vh^ = 

£ Jo Jq 



± 


- / ro,s(piO) - 


In 


fl 

Jo Jn 


■dt(p, 


which obviously converges in the limit e ^ 0. 

In order to deal with the transport and the capillarity terms, we want to use the structure of 
the system. Therefore, first of all we need to introduce a regularization of our solutions. 


4.3.1 Regularization and description of the oscillations 

Let us set 14 := Pe^e- We can write system ([12]) in the form 


(37) 


sdtVe + div 14 = 0 

edtVe + (^c{x^)e^ x ^ + V(ld - A)r,'^ 


e/e 


19 



where we have defined by the formula 

(38) /e := -div {peUe®Ue) + V d\\ {peDUe) - 

- ^v(n(/,,) - n(i) - n'(i) {p, - 1 )) +^{Pe- 1) VAp,. 

Equations (I37p have to be read in the weak sense, of course. In particular, from writing 

(/e,'^) := jJ^PeUe ® Ue : V ijj - V PeDUe : V ijj - ^/\peV Pe ■ iIj - 

1) APedivV’ + ^{jl{pe) - n(l) -n'(l) {pe - l))divV>^dx 

fe ■ + fe : VV' + /I • V’ + fe div V' + fe div V') dx 

and by a systematic use of uniform bounds (recall Propositions 14.1] and HiBI) . we can easily see 
that (/e)^ and (/|)^ are uniformly bounded in L“(L^), and so is (/|)^ in L^{L^); finally, (/|)^ 
and {f^)^ are bounded in L|n(L^). 

Therefore, we deduce the uniform boundedness of (/e)^ in the space + iy“^’^(n)), 

and then in particular in Ly(if“®(0)) for any s > 5/2. 

Now, for any M > 0, let us consider the low-frequency cut-off operator Sm of a Littlewood- 
Paley decomposition, as introduced in (15011 below, and let us define 



re,M := SmTc and := S'mK . 


The following result hods true. 

Proposition 4.6. For any fixed time T > 0 and compact set K d FL, the following convergence 
properties hold, in the limit for M —)• -|-oo; 


(39) 


SUP£>o Ike '^e,M\\L'^{H={K))r\L‘^{m+‘‘{K)) Vs<l 

SUP£>oll'f^£ - ^0 Vs>0. 

Moreover, for any M > 0, the couple {fie^M j ye,M) satisfies the approximate wave equations 


(40) 

where {fe,M 


{ £ dtrs,M + div Ve,M = 0 

edtVe,M + {c{x^)e^ x + V(ld - A)r£,M) = £ fe,M 
)^ and [ge,M)^ families of smooth functions satisfying 


+ ffe.M , 


(41) 


^^Pe>o\\fs,M\\L'^[H^{K)) — Vs>0 

SUP £>0 ||5e,Mlk2^(j:^i(^)) —>0 for M-Foo 

where the constant C{s,M) depends on the fixed values of s >0, M > 0. 


Proof. Keeping in mind the characterization of spaces in terms of Littlewood-Paley decompo¬ 
sition (see Appendix lA.ip . properties (I39D are straightforward consequences of the uniform bounds 
established in Subsection KT\ 

Next, applying operator Sm to dSZl) immediately gives us system (sni), where, denoting by 
[P, Q] the commutator between two operators V and Q, we have set 

fe,M ■= SMfe and ge,M ■= [ c ( x ''), 5 m ] (e^ X Ve) . 
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By these definitions and the uniform bounds on (/g) , it is easy to verify the first property in 
(HU). As for the second one, we need to generalize the argument of Lemma 3.3 in [25], because 
here we have less available controls on the family (14) . 

First of all, by uniform bounds and Lemma lA.7l we get 

sup||5(£,M|li,2(i2) < C2-’^ . 

e>0 ' 

As for the gradient, for any 1 < J < 3 we can write 

djQsM = [c, SM]dj (e^ x I4) + [djC, Sm] (e^ x I4) . 

In order to control the former term, we use Lemma lA.81 with p 2 = q = 2 and pi = 1. Recalling 
that, by (l26]) . (D14)^ C Ly(L)^^), for any compact K C 0, we get 

sup II [c, SM]dj (e^ X K)||i 2 ^(i 2 (^)) < C2-^ . 

For the latter term, instead. Lemma lA.91 gives us 

sup II[«9jC, Sm] (e^ X 14)11 2 /^ 2 ^ < Cfj,{2~^). 

£>0 ^ 

In the end, choosing rj{M) = max |2“^,/r(2“^)| (which goes to 0 when M —)• +00), we get 

sup ||5e,M||£,2 ^ < Crj{M) 

for a suitable constant (7 > 0, and this completes the proof of the proposition. □ 

We also have an important decomposition for the approximated velocity fields. 

Proposition 4.7. The following decompositions hold true: 

14,M = and DVe,M = ©e.M + 

where, for any compact set K <zTt and any s > 0 one has 

f + ll®e-^llL2,(L2(i<')) “ 

|||V£,m||^^(^.(^)) + ll^e,M||^2(^.(;^)) < C{K,S,M), 

for suitable positive constants C{K), C{K, s, M) depending just on the quantities in the brackets. 

Proof. By definitions, we immediately have 

Ve,M = SMiPeUe) = Sm{pI^'^ U^) - S m Ps Ue) ■ 

Thanks to the uniform bounds established in Subsection im the first decomposition and the 
related estimates are easy to be verified. 

Let us now take a space derivative of 14,Mj splitted in accordance with the previous identity. 
Thanks to spectral localization, the second term do not present any problem: indeed, for any 
1 < j < 3 one has 

dj ((“ 1 ) Pe Ue) = dj - - Pe Ue + - - dj [pe uf) 

= VP^ Ue + - - dj {pe Ue) . 

Zi s. s 
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For the former term, one uses Propositionand the fact that (by Sobolev embeddings) G 

L|,(L®); for the latter, instead, one has to take advantage of the estimate \^/p—l\ < |/9—together 
with property (I26p . 

As for Sm Ue )) instead, we have to proceed carefully. More precisely, for any 1 < J < 3 
we start by writing 

djSu Ue) = ^Sm {yfFe Ue djPe) + Sm [pI^"^ djUe) 

3 

= -^SMi^/fhUedjPs) + SMiVP^djUs) + S'm ((Ps - 1) • 

Recalling now that (^y/p^Due)^ C L|,(L^(n)), we can set = SM{^Jlh^jUs)■ Indeed, since 

(e“^ {Pf. — 1))^ C L^(L®(n)) and C {^Lp‘ {QPj), the other two terms are of order 

e. So we can include them into the remainder V^^m- 

Hence, the proposition is now proved. □ 


4.3.2 The capillarity term 

First of all, let us deal with the surface tension term in (I36F Notice that it can be rewritten as 
fo ■ '0 , for any smooth test function •0. 

Thanks to the next lemma, we reconduct ourselves to study the convergence in the case of 
regular density functions. 

Lemma 4.8. For any ip G iD([0, T[ xll; , we have 


lim lim sup 

M^+oo e^o 


n Are Vre ■ Ip dxdt - / Ar^^M ■ ipdxdt 

! Jo Jn 


= 0 , 


Proof. Let us fix M > 0: we can write 



Are Vre - Ip = [ [ Are Vre,M ■ ip + [ [ Ar^ V(id - SM)re ■ V’ ■ 

/o Jn Jo Jn Jo Jn 

By the uniform L|,(L^) bounds on the family (Are)^ and the first property in (13911 . we get that 

rT 


0 JQ 


Are V(ld - SM)re ■ ip 


< C6{M) 


for some positive constant C, independent of e and M, and some function 5{M) such that 
6 {M) —> 0 for M —?> +00. 

On the other hand, the former term in the right-hand side of the previous identity can be 
written as 



Are Vre,Ar • ip 


0 Jn 


Are,M Vre,M • '0 + / / A(ld - SM)re ■ ip ■ 


0 Jn 


For the last term, using that the operator Id — Sm is self-adjoint, we can estimate 


[ [ A (Id - S'M)re Vre, AT •-0 
Jo Jn 


< C IIAre 




I (Id - SM){yrs^M ■ 


L2(L2) 


and this quantity converges to 0 for M —)> -|-oo, thanks to Lemma [A.31 point [in). Indeed, it is 
enough to notice that (f^iyr^^M • ^))e is bounded in L|.(L^), uniformly both in e and M. □ 
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Thanks to Lemma l4.8[ for any ijj G P([0, r[ xO; H Ker^ we have to consider the con¬ 
vergence of the term (pay attention to the signs) 

n Ar^^M■ ipdxdt = - (id - A) • V'dx dt 

! Jo Jn 


+ / / -ipdxdt. 

Jo Jn 

Notice that = V 1‘J'- therefore, since ^ is divergence-free, by integration by 

parts we get that the latter item on the right-hand side is identically 0. 

Therefore, in the end we have to deal only with the remainder 

(42) - [ [ (Id - A)r£,M Vre,M •'0 = ~ [ [ - Ah) {re,M)'^h{re,M) ■ - 


0 Jn 


0 Jn 


((Id - A)r^^M Vr^^M) • V', 

0 Jn 

where we used the notations introduced in ([2]), setting the mean-free part of r^^M- 

4.3.3 The convective term 

In the present paragraph we will deal with the convective term. Once again, the first step is to 
reduce the study to the case of smooth vector fields 14,M- 


Lemma 4.9. For any if; G iD([0, T[ xll; M^), we have 


lim lim sup 

M^+oo £—>.0 


a p^U^ ® ■.V'll) dx dt — / / Ve^M ^Ve,M dx dt 

I Jo Jn ’ ’ 


= 0 . 


Proof. First of all, we can write 

n PeUe (g) Ue : V'f/’ = 

! 



0 Jn 


f \fPe ~ 1 

PeUe® y/PeUe -.Vlf - e j / - PeU^ ® ■. V'if 

Jo Jn ^ 


Notice that the latter integral in the right-hand side is of order e: this is a consequence of the 
uniform bounds (^^/pfu^)^ C L^^L"^) and {^/pf — 1))^ C L^{L°°) for any p G [2,4[ (recall 
Proposition 14.31) . The former one, instead, can be decomposed again into 

n peUe® yfplUe -.Vlf = ® y/fk Ue ■ + / / {U - SM)Ve ® ^/pf Ue : V'lf 

! Jo Jn Jo Jn 

= [ [ Ve,M®Ve:Vlf + 0{e) + [ [ {Id -SM)Ve®^eUe:Vf;, 

Jo Jn Jo Jn 

where, in passing from the first to the second equality, we have argued exactly as before. 

Let us focus on the high frequency term first: we can write 


f f {Id - SM)Ve® ^/pfUe :Vif = -e f [— - - {id - SM)Ve ® ^/fkUe : Vlf + 

Jo Jn Jo Jn £ 

+ [ [ {Id -SM)Ve®pl/^Ue:Vlf. 

Jo Jn 

Notice that the former term is 0(e) (by uniform bounds again). On the other hand, (I27D tells us 
that {pe^'^Ue'Vif) ^ is uniformly bounded in PiKPia)): then, since (V^)^ C L‘^{L‘^) and the 
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operator Id — Sm is self-adjoint, by Lemma [A.31 we gather that the latter one is arbitrarly small 
for M large enough, uniformly in e. In the end, we obtain that 


lim lim sup 

M^+oo 



0 Jn 


(Id - 5 m)I4 ®y/lhUe : VV' 


= 0 . 


It remains us to consider the integral 


IL I4,m ® K,m : VV' -h [ j I4,m ® (W - 5'm)I4 : VV' 


0 JO, 


0 Jo. 


and prove that the last term in the right-hand side of the previous relation is small. We notice 
that, by the decomposition of Proposition 14.71 we have 

||D(I4,m : <Ci+ eC2 iM), 

for some constant Ci > 0 independent of e and M, and C 2 {M) just depending on M. Therefore, 
using again the simmetry of the operator Id — Sm and Lemma IA.3I it immediately follows that 


lim lim sup 

M^+oo 


r [ K,M®(ld -5 m)K:VV' 
Jo Jn 


= 0 . 


This concludes the proof of the lemma. 


□ 

Now, recall equation (I36p : paying attention once again to the right signs, by the previous 
lemma we have just to pass to the limit in the term 


/ [ VeM ® K,m : VV' = / / div (14 ,m ® 14,m) • 
JO Jn Jo Jn 


f /div4( 

Jo Jn ^ 


^£,m) ® (K,Af) )•'*/’ + 


f-T 

lo Jn 




-0 



0 Jn 


+ 'T'em) ■ ■ 


For notational convenience, from now on we will generically denote by TZs^m any remainder, 
i.e. any term satisfying the property 


(43) 


lim lim sup 

M^+oo £—>0 


f! 

Jo Jn 


^£,M • ijjdxdt 


= 0 


for all test functions ijj £ 'D([0, T[ xll; M^) n Ker^. 

Handling T^M‘ Since we are dealing with smooth functions, we can integrate by parts: we get 

T^m = div;,((y^^)®(y^^)) = div,((y,"^)) 


= div,((y^Aj)) {V^,m) + 2^" V J ^ KmV- 

Notice that we can forget about the second term: it is a perfect gradient, and then it vanishes 
when tested against a function in the kernel of the singular perturbation operator. 

For the first term, we take advantage of system (1401) : averaging the first equation with respect 
to and multiplying it by we arrive at 


diVh((14^M)) (K^m) — ~ ^ dt{re,M) — ^£,M + £ (^e,M) 5* (I4 ^m) ) 
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since £ dt[{rg^M) ^ remainder in the sense specified by relation (1431) above. We use now 

the horizontal part of (1401) (again, after taking the vertical average), multiplied by paying 

attention to the signs, we get 

e {re,M) “ (^^m) - ^h){reM) + ^ U^m) + (^m) 

= - c(x^) {re,M) + (M - ^h) {Tsm) ^hircM) + , 

where we used also the properties proved in Proposition 14.61 and we included in the remainder 
term also the perfect gradient. Inserting this relation into the expression for 7^^^, we find 

(44) %]m = (curl,(y^M) - c(x") {Ve%)^ + (H - A,) {r,J V,(r,^) + 7^,,A^ . 

In order to deal with the first term in the right-hand side, the idea is to decompose in the 
orthonormal basis (up to normalization) {V/iC, V^c}. Of course, this can be done in the region 
when V/iC is far from 0: therefore, we proceed as follows. 

First of all, let us come back to system (jlHl) for a while: we take again the vertical average of 
the equations and we compute the curl of the horizontal part, finding 

edtcm\h{V,^^M) + div/^ (^c(x^) = ecuilhif^^M) + curlh{g^M)- 

On the other hand, from the first equation multiplied by c, we get 

edt (c(x^) (re,M)) + div,, (^c(x^) (V;^m)) = (K^) ' , 

and this relation, together with the previous one, finally gives us 

(45) edt (curlhiV^^M) - c(x^) (r^^M)) = ecurl/,(/^'(^) -h curlh{g^M) + (K^m) ' Vfec(x^). 

Notice that, thanks to Proposition 14.61 there exists a function r/ > 0, with g{M) —)• 0 for 
M 4-00, such that, for any compact iP C fl. 


(46) 


sup 

£>0 


cmlhig^^M) 




< v{M ). 


Then, fixed a 6 S C((°(M^), with 0 < b{x^) < 1, such that 6 = 1 on {|x^l < l} and 6 = 0 on 
{\x^\ > 2}, we define 

bMix’^) := b(^{g{M))-^^^Vhc{x^)) . 

Now we are ready to deal with the first term in the right-hand side of (I44p . For notational con¬ 
venience, we set ■= curl/i(I4^^) — c(x^) {ts^m)- On the one hand, using the decomposition 

and the bounds established in Proposition 14.71 we deduce that, for any compact iP C fl. 


bM Xe,M 


Li([0,r]xX) 


< eC{M) + C 

< £C{M) + C {c\x'^ Vhc{x^) 


< 2 


1/6 


Therefore, thanks to hypothesis (1171) . we infer that this term is a remainder, in the sense specified 
by relation ()43ll . On the other hand, for V/jC far from 0, we can write 


{l-bM) 


(1 — bu) ^e,M 


'{Vs 


£^m) 




|Vhc( 


V^c + 
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= (1 - ^m) 


e,M 


■V^c + —-V/^c). 


iVhcr 


iv,,cr 


We observe that the latter term in the right-hand side is identically 0 when tested against a 
ip G Ker^ (because div/i'0^ = div/j (c?/;^) = 0, and then ip^ ■ V/^c = 0). For the former term, 
instead, we use the expression found in (1451) : we get 


(1 - 6m) 


M 


(^^M)-VhC 


V^C = 


|V/^c|2 

(1 - 6m) d4,M 

iv,ci^ --iv,cr 

e (1 — 6 m) i9t |d4,Mp „_L (1 — &m) d4,M 

— - 


<.Fiy V7-Lr (l-^M)d:'£,M 
^ Ot^e,M V/j C — 


2 |V;,cr 


iVhd 


ecurlfe(/^^Af) + curl/,(5^^j\^)) V^c 
ecni\h{fe,M) + curU(5^_^)) V^c, 


which is again a remainder TZs,M^ thanks to Proposition 14.71 and property (|46p . 

In the end, putting all these facts together, we have proved that (paying attention again to 
the right signs) 

(47) TIm = (Id - A;,) V,(r,^) + 77 ,,m • 

Notice that the first term in the right-hand side exactly cancels out with the first one of (jl2]). 
Dealing with TZ^m* consider the term 77 ^m- exactly as done above, we can write 


■trh 

^e,M 


) + {cmlhV^j^ (K^m) ) 


T!,m = div, ® = (div^V^M) Vsm) + 2 

Let us focus on the last term for a while: with the notations introduced in Q, we have 


X 


curlF,,M = 


M 


and 


where we have defined := (v^m^ ~ ^3 1^°'- these quantities, from the momen¬ 

tum equation in (14011 (where we take the mean-free part and the curl), we deduce 


curll4,M ) = curhV;,M = ^e,M > 


(48) 


- cKV = (<93 ^ curl (e/,,M + 5e,M 


'^Im + divfc(cy^^) = curU + 9^ 




Making use of the relations above and of Propositions 14.61 and 14.71 we get 

X 


curlfe^'^M (I^V 


rh 


= OJl 


'e,M \^e,M 


= -dt 


't I 


7)3 I n3 

^e,M ~ ~‘^£,M 


^3 ^ curl ( e fe,M + ge, 


M 


= - - w 


'e,M 


dti^, 


e,M 


77, M = 




(m) div/i(cl4^M) + T^£,m ■ 


Hence, including also the gradient term into the remainders, we arrive at the equality 

= (div,(H^^) (v,% + (w^,mY)) + (i (w^,mY I?m-V,c) + 77,,M, 
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which can be finally rewritten in the following way: 


T; 


e,M — (div )) “ 


“ ^e^M'^/iC) + 7le,M ■ 

The second term on the right-hand side is actually another remainder. As a matter of fact, direct 
computations yield 


( K,M + [^e,M j 


rh 


= K 


'31 '^£,Af yje,M 
1 


KV + w, 


' e,M 


-VlMdAV,^,M+{W, 


'e,M 


= T^e,M “ 2 


w3 

^e,M 


= T^e,M ■ 

As for the first term, instead, we use the equation for the density in (1401) to obtain 
div Ve,M + (^e^m) ^ = -^ dtr^^M i^e,M + {^e,M 

= Tle,M + ^’’e.M 

Now, by equations (|48p and (I40p again, it is easy to see that 




■M 


= '^e,m - r£,M V(ld - A)re,M = '^e,m + Vr^^M (id - A)re,M , 


and therefore we find (with attention to the right sign) 

(Id - {w^mY • V,c) + 7^,,M. 

Now we have to deal with the second term in this last identity. Once again, we take advantage 
of the decomposition along the basis {V/jC, V^c}. For simplicity of exposition, we omit the 
cut-off away from the region {V/^c = 0} and we just give a sketch of the argument, since it is 
analogous to what done above for 
First of all, we write 


Km 


Ve.M-^hC = -KM-^hc]{W,%-VhC 


iVhCp 


+ W, 


' e,M 


■ VftC 


VfeC 

iv,.cr 


As before, the last term in the right-hand side vanishes when tested against a smooth ijj E Ker^. 

T,m 


Next, we obtain informations on V!}j^ ■ V/jC from the first equation in 




1 


]m ■ VftC = ^ [edtW^^M - (^^3 ^ curl [e fe^M + Qe, 


■M 


VftC. 


Therefore, we obtain that 

1 
c 


Km 




2H 


dt 


w!^ 


e,M 


VfeC 


|Vhc|2 


dn ^ curl ( e fe,M + g£,M 


VhC 




IVftcl 


2 ’ 


which is obviously a remainder in the sense of relation 

In the end, we have discovered that (paying attention to the right sign) 

(49) Tf j^ = (Vre^M (Id — A)re,M) + ■ 

Notice that the density-dependent term exactly cancels out with the latter item in 
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4.3.4 The limit equation 


Let us sum up what we have just proved. In order to pass to the limit in equation (|36p . we needed 
to treat the non-linearities coming from the capillarity term and the convection term. 

For the former, after applying Lemma 14.81 we have arrived at relation 


- [ [ (Id - A)re,M Vre,M • "0 = ~ [ [ (Id - A/i)(re,M) V/j(re,M) •-0 - 
Jo Jo. Jo Jn 

- f [ {{Id - A)re,M Vre,M) • Ip , 

Jo Jn 

On the other hand, for the latter we exploited Lemma [4.91 and, after manipulations, we have fonnd 

~ [ [ ® 14,M : VV’ = / [ {T}^m + 'T'sm) ■ V' 

Jo Jn Jo Jn 

(^(Id - Ah){re^)Vh{rsj^) -|- (Vfe,M (id - A)re,M) + • -0, 


cT 

lo Jn 


where we used also relations (I47|) and (|49p . Therefore, thanks to the special cancellations involving 
the density-dependent terms, in the end we have proved the relation 

[ [ (Id - A) r£,M Vr£,M • '0 - 14,M <8) 14,M : • 0 , 

Jo Jn\ / Jo Jn 

which immediately implies, together with Lemmas 14.81 and 14.91 that 


lim lim 


1 


Ape Vpe -0 — Pe Ug O Ug : V0 \dxdt = 0 . 


M —>-+00 £—^-O Jq J\^S 

Then, thanks to the previons compntations, we can pass to the limit in the weak formnlation 
of onr system: we obtain 


>0 Jn 



—u • (9t0 — r(9t(ld — A/i)0 -|-: V0^ dxdt = J (uq • 0(0)-|-ro(ld — A/j)0(O))(ix 


for any (0, 'll;) test fnnctions belonging to the kernel of the singular perturbation operator A. 
Recall that, in particular, this implies the relation c0^ = V^(ld — Ah)cj). 

Fnrthermore, we recall that also (r, tt) € Ker.A: then divu = 0, u = (u^,0) and cu^ = 
V^(ld-A,)r. 

Setting X{r) = (id — A^jr and 0 = (id — A/i)0, and using that all the fnnctions depend 
jnst on the horizontal variables, straightforward compntations yield to 


u ■ dt'ip dx dt = 


lo 

cT 


div/j( -^VhX{r) ] dtcpdx^dt 


[ f Du'.Vipdxdt = / f iDc(X(r)) : Vh{c~^ dx^ dt 

Jo Jn Jo Jr'^ 

= zz / f Tt,:(A(r)) : 'Sc{(p)dx^dt = i/ / f o 'S^{X{r)^ cpdx^ dt. 

Jo ./r 2 Jo Jr'^ 

Inserting these equalities into the previous relation completes the proof of Theorem 12.81 
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4.4 Final remarks 


We conclude by making a few comments about our hypotheses on the regularity of the rotation 
coefficient, namely c G with V^c G 


• At the level of uniform bounds (see Section [3]), we asked for V/jC G L°° in order to close the 
estimates for the BD-entropy (recall Lemma 13.61) . both in the case of vanishing and constant 
capillarity. Nonetheless notice that, under an additional L|,(L^) bounds for the family of 
velocity fields (ue)^, the control of the rotation term in (1191) is straightforward for c G L°°: 


1 

£ 



X Us ■ V ps dx dr 


< CVt ||c||io 


l|Vp, 






where the right-hand side is uniformly bounded by classical energy estimates (see Corollary 
ESI). Remark that the property (us)^ C Ly(L^) can be deduced, for instance, when an 
additional laminar friction term is added to the momentum equation (see e.g. El). 

• The W^’°° regularity seems to be necessary in the compensated compactness argument: in 
particular, we exploited it for resorting to the equation for the vorticity and for decomposing 
some vector fields along {V/jC, V^c}. 

Notice that, by Proposition 14.41 V^c is parallel to V/j(Id — A/i)r; on the other hand, we 
have no available informations for V/j(Id — A/i)r: in particular, it is not clear for us how to 
avoid the non-degeneracy condition (HZD. 


• Finally, the requirement V/jC G (for some admissible modulus of continuity p) is funda¬ 
mental in order to get (|46l) . which is a key property in our proof. Notice that paraproduct 
decomposition does not seem to help in gaining anything at this level. This is the main 
reason why imposing some regularity on seems to be necessary for proving the result, 
and it seems not to be possible to go below the Lipschitz threshold c G 


Related to the non-degeneracy condition, let us remark also that there is a huge gap between 
the case of constant rotation axis (for which Vc = 0 everywhere) and the case of variable rotation 
axis considered here, for which we need hypothesis m- It seems very likely that this difference 
is just “artificial”, and it depends on the techniques used in the proof. It could be interesting to 
find a different approach to the problem, in order to being able to treat coefficients which can be 
constant (i.e. whose gradient can vanish) in a region of non-zero Lebesgue measure. 


A Appendix 

We collect in this appendix some well-known results about Littlewood-Paley theory and admissible 
moduli of continuity, and the proof of some technical lemmas. 

A.l Fourier Analysis toolbox 

Let us recall some tools from Fourier Analysis, which we exploited in our proof. We refer e.g. to 
[2] (Chapter 2) and |27| (Chapters 4 and 5) for the details. 

For simplicity of exposition, let us deal with the case; however, the construction can be 
adapted to the d-dimensional torus T'^, and then also to the case of x T'^^. 

First of all, let us introduce the Littlewood-Paley decomposition, based on a non-homogeneous 
dyadic partition of unity in the Phase Space. 

We fix a smooth radial function y supported in the ball i?(0, 2), equal to 1 in a neighborhood 
of i?(0,1) and such that r i—>■ x(r e) is nonincreasing over M_|_ for all unitary vectors e G Set 
p{0 = X (0 - X (20 and := for all j > 0. 
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Aj := (p{2-W) if j > 0. 


The dyadic blocks (Aj)jeZ ai'e defined b}@ 

Aj := 0 if j < —2, A_i := x{D) and 
We also introduce the low frequency cut-off operators: for any j > 0, 

(50) SjU := AkU. 

k<j-l 

The following classical property holds true: for any u G S', one has the equality u = Ylj 
in the sense of S'. We will freely use this fact in the sequel. 

Spectrally localized functions have nice properties with respect to the action of derivatives. 
This is explained by the so-called Bernstein’s inequalities. 

Lemma A.l. Let 0 < r < R. A constant C exists so that, for any nonnegative integer k, any 
couple {p, q) in [1, -|-oo]^ with p < q and any function u G L^, we have, for all A > 0, 

suppuC B{0,XR) ||V^u||l, < ||rx||LP ; 

suppu C U G I rA < Id < RX} C-^~^ ^’"\\u\\lp < \\^’"u\\lp < X^\\u\\lp • 


By use of Littlewood-Paley decomposition, we can define the class of Besov spaces. 

Definition A.2. Let s G M and 1 < p,r < -|-oo. The non-homogeneous Besov space is 
defined as the subset of tempered distributions u for which 


•= 




feN 




< - 1-00 . 


Besov spaces are interpolation spaces between the Sobolev ones. In fact, for any A: G N and 
p G [l,-|-oo] we have the chain of continuous embeddings B^i ^ ^ ^poo- Moreover, for 

all s G M we have 2 = , with equivalence of the respective norms. 

We report here some statements which we used in our analysis: their proof are given in m- 
We refer again to [2] and m for more details and more general results. 

Lemma A.3. (i) For 1 <p <2, one has ||/||l 2 < C'(||/||i,p ||V/||i 2 ). 

(a) For any 0 < 5 < 1/2 and any 1 < p < -|-oo, one has 


II/IIl- < c[\\fUp + llv/lli-V')-'' llvVIliV'^^'). 


(in) Let 1 < p < 2 such that 1/p < 1/d + 1/2. For any j G N, there exists a constant Cj, 
depending just on j, d and p, and going to 0 for j —)■ -|-oo, such that 

II(Id - 5,)/IL2 < C, ||V/||^o^ . 

In particular, ifS/f = V/i -|- V/ 2 , with V/i G 00 V /2 G Bp^^, then 

||(Id - Sj)f\y < Cj (liv/illgo^ + llVMIso^) , 

for a new constant Cj still going to 0 for j —)■ -|-oo. 

^Throughout we agree that /(D) stands for the pseudo-differential operator u 1—^ Tu). 
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Finally, let us recall some notions of homogeneous dyadic decomposition. Namely, one can 
rather work with homogeneous dyadic blocks {Aj)j^z, defined as 

Aj := D) for all j £ Z. 

Then, we can introduce the homogeneous Besov spaces Bp j. by the property 




£r 


< +00 


(plus some conditions on low frequencies). We do not enter into the details, for which we refer 
to Chapter 2 of [2]. Let us however recall rehned embeddings of homogeneous Besov spaces into 
Lebesgue spaces (see Theorem 2.40 of [2]). 

Proposition A.4. For any p £ [2,+oo], one has the embeddings Bp 2 ^ and LP' B^, 2 , 
where p' is defined by the eondition 1/p' = 1 — 1/p. 


A.2 Admissible moduli of continuity 

In this paragraph we recall some basic definitions and properties about general moduli of conti¬ 
nuity. We refer to Section 2.11 of [2] for a more indeep discussion. 

Definition A.5. A modulus of continuity is a continuous non-decreasing function p : [0,1] —)■ 
M+ such that p{0) =0. 

It is said to be admissible if the function T^, dehned by the relation 

Tp{s) := sp{l/s ), 


is non-decreasing on [1, -|-oo[ and it verihes, for some constant C > 0 and any s > 1, 


r+oo 

/ a~'^Tp{a)da < Cs~^Tp{s). 
J S 


Given a modulus of continuity /r, we can define the space C^(M“) as the set of real-valued 
functions a £ L°°(M'^) such that 

I I \a{x + y) - a{x)\ 

|a|c„ := sup -TT-r^- < -|-oo . 


h|G]o,i] 


Fi\y\) 


We also define ||o||c^ := ||a||L°° + |o|c^- 

On the other hand, for an increasing T on [1, -|-oo[, we define the space By{M.'^) as the set of 
real-valued functions a £ L°°(M'^) such that 


I 11 ^ ® 11 

a\Br ■= sup———— < -hoo, 
j>0 1 ) 


where Sj is the low-frequency cut-off operator of a Littlewood-Paley decomposition, as introduced 
above. We also set ||a||Bp := ||a||L°° + 

One has the following result (see Proposition 2.111 of [2]). 

Proposition A.6. Let p be an admissible modulus of continuity. Then = Br^(Sfi), and 

the respective norms are equivalent. Moreover, for any a £ C^(M'^) one has 

W^jaW^oo < Cp{2-^) 

for all j > —1, where the constant C just depend on ||a||c^. 
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Now we want to present a commutator lemma, which is fundamental in the proof of Proposition 
14.61 especially for property (1411) . 

First of all, let us recall the classical commutator estimates (see Lemma 2.97 of El). 

Lemma A.7. Let 6 £ such that (l + | • |)0 € L^. There exists a constant C such that, 

for any Lipschitz function I G and any f G LP{W^) and for all A > 0, one has 

\\[9{X-^D),£]f\\^^ < CX-^ \m\L^ II/IUp. 

Going along the lines of the proof, it is easy to see that the constant C depends just on the 
norm of the function |x| k{x), where k = denotes the inverse Fourier transform of 9. 

Let us give a slight variation of the previous lemma. For simplicity, we restrict our attention 
to the case of 9 in the Schwartz class this will be enough for our aims. 

Lemma A.8. Let 9 € and {pi,P 2 ,q) G [Ij+oo]^ such that l/q = 1 + l/p 2 — 1/Pi- Then 

there exists a constant C such that, for any f G any I G and all A > 0, 

\\[9{X-^D),e]f\\^,^ < CX-^ \m\L^ II/IUpi. 

The constant C just depends on the L^ norm of the function |x| k{x), where k = J-^^9. 

The proof follows the arguments used for the classical statement, with no special novelties. 
Hence, we omit it. 

Let us consider now less regular functions i. 

Lemma A.9. Let 9 G C^(R'^) be as in Lemma IH.71 and let p, be an admissible modulus of 
continuity. Then, there exists a constant C such that, for any function I G (^^(R'^) and any 
f G LP(R'^) and for all X > 1, one has 

\\[9{X-^D),i]f\\^^ < C/r(A-i) \i\c^ ll/lliP. 

The constant C only depends on the L^ norms of the functions k{x) and |x| k{x). 

Proof. As in the proof of the classical result (see Lemma [A. 71) . we can write 

[9{X~'^D),i] f = X‘^ [ k(X(x - y)) f{y) {£{x) - £{y)) dy . 

Jm'i 

Remark that the previous integral is actually taken over R'^ \ {x}, so that we can multiply 
and divide by y(\x — y\). Making the seminorm appear, thanks to Young inequality we are 
reduced to estimate the quantity 

ALltOOMd-Dlli. = A" / mxz)^{\7.\)d7.. 

Let US split the previous integral according to the decomposition R'^ = {|2:| < A“^}u{l 2 ;l > A“^}. 
For the former term, since y is increasing we have 

X"^ [ \k\{X z) fj.{\z\) dz < fj.{X~^)\\k\\Li . 

i|z|<A-i 

For the latter term, instead, we make the non-decreasing function F^ appear, and we estimate 
A'^ [ \k\{X z) y{\z\) dz = X^ j \k\{Xz)T ^{\z\~^)\z\dz 

< cr^(A)A-i II I-IMOIIli < Cy{X-^). 

The lemma is hence proved. □ 

Obviously, an extension of the previous result, in the same spirit of Lemma lA.81 holds true. 
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A.3 On the BD entropy strnctnre 

We give here the details of the proofs of some technical lemmas about BD entropy estimates. We 
start by proving Lemma [3.41 

Proof of Lemma \3.4\ First of all, by Lemma 2 of [7] we can write 

(51) - -f- [ /? |V log + [ Vdivtt • V/9 + [ pDu : V log p (g) V log p = 0 . 

2 dt Jq Jq 

Next, we multiply the momentum equation by vSIpjp and we integrate over D: we hnd 


(52) V / (ptu + u ■ Vu) • Vp + / Du : f V^p 

Jfi Jo. V 


--Vp® Vp ) + 
p 


+ — / X tt • Vp + 
Ro 




Now we add m, multiplied by to f|52p : we end up with 

,2 


+ Wei 




= — u 


' / dtu ■ p — / Vdivri • Vp — {u ■ Vu) ■ Vp — / Du : V^p. 

Jq Jq Jq. Jq. 

From this relation, thanks to the mass equation and the identities 

— / « • Vdiv (pti) — / {u-Vu)-Vp = / pVu:^Vu 

J Q, J Q. J Q, 


— / Vdivrt-Vp— / Du-.V^p = 0, 

Jvt Jq 


we deduce the equality 

/JP'(P) + Kip)) IV^P <ix + 


We 


in 


Fr 


In 


+ — [ X u ■ Vpdx = — V ^ [ u ■ Vpdx + z/ / p Vu : ^Vudx . 

Ro Jq dt Jq Jfi 

Notice that we can rewrite this last relation in the following way: 

^4: [ P \u + vVlogpl^ dx + [ IV^pl^ dx + ^ X u -Vpdx + 

2 dt Jq Vie Jn Ro Jn 

+ ^ [ (P'iP) + PciP)) |VVp|^ dx = ^4 / p\u\‘^dx + u [ pVu : ^Vudx. 

Then, we conclude by integrating with respect to time and using Proposition 13.21 □ 

Now, let us switch our attention to the proof of Lemma 13.51 
Proof of Lemma \3 . 51 By Lemma [A.31 item (i), we easily deduce 

\\P-M\l'^{D) < C (^\\p - + ( 1 - 12 ( 7 )) II^PllLf=(L 2 )^ , 

and the first estimate immediately follows from Proposition 13.21 and Corollary 13.31 

Let us now focus on the second estimate. By Lemma lA.31 item (ii), for any 0 < h < 1/2 and 
any 1 < p < +00 we can write 


(53) 


||p-i|lip(^^) < Cp{\\p-i\\Lri.D)t + dT^ . 
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Now, since p < 4/(1 + 26), we can apply Holder inequality to make the Lj(L^) norm of 
appear. More precisely, using again the bounds of Proposition 13.21 and Corollary 13.31 we have 

< (We)P((V2)-^) ^1-1/. (||vV||'2(i.))'^' , 

where we have set q := 4/((l + 26)p). Putting this inequality into (I53p and using also the bound 
of the first part, we finally get the desired estimate. □ 


Let us conclude this section by showing the proof of Lemma 13.61 Notice that we are not able 
to exploit the presence of the cold pressure term at this level. 

Proof of Lemma \S.6l The first inequality is trivial: we have just to write x u ■ S/p = x 
{^/pu) ■ Vy^, and then apply Holder inequality to the integral over H and Proposition 13.21 
Let us then focus on the estimate in (ii). First of all, we write 

(54) f f ce^xu-Vp= f f ce^ X yPpu ■ Vp + f f 

Jo Jn Jo Jo, Jo Jn 

Now we perform an integration by parts in the latter term: denoting by cj = V x u the vorticity 
of the fluid, we get 

[ [ ce^ X yf^u-Vp(^-l] =11 cpuj^{^-l) + 

JoJn \y/P / JoJn 

+ [ [ pu- V^c ( Vp - 1) - f [ ^ X • Vp xu-Vp, 

JoJn JoJn^ \y/P / JoJn^ 

which in turn implies 

//ye“x^„.Vp(T-l) = IIJ^c^UV-p-l) + 

+ ^/ / ce^X'«-Vp+^ / [ pu- {y/p-l) . 

^ Jo Jn -J Jo Jn 

Combining now (fMl) with this last relation gives us 


ce^ X yJpu ■ Vp ( - 1 

\Vp 




0 Jn 


Ro 


ce^xu-Vp = — / cpoj^{^/p-l) + 


0 Jn 


Let us start by considering the second term on the right-hand side: on the one hand, by use 
of Holder inequality and Proposition 13.21 one gets 


( 56 ) 


u 

Ro 



X ^/pu ■ Vp 


< 



||Vp||j ;^2 dr < Cut 


VWe 

Ro 


On the other hand, if 7 = 2, by use of Young inequality we can also write 


( 57 ) 


u 

Ro 



X ^/pu ■ Vp 
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For the vorticity term, we start by observing that \y/p— Ij < |/9 — 1|: therefore, thanks also to 
Corollary 13.31 and Lemma [3.51 we can write the estimate 

ct /* 

CpU^ {y/p-l) 


Ro 


0 Jo, 


V 

Ro 

Cv 

Ro 


^ ^ C \Wp\\l‘^(l°°) ’ 


where, for notation convenience, we have set 

C = C(Fr, We) := (ft + (l - 12 ( 7 )) VWe) . 

In order to control the last factor in the right-hand side, we take advantage as usual of the 
decomposition = 1 -|- (y^ — 1); then, applying Lemma 13.51 with e.g. 5=1/4 and p = 2 (and 
so g = 4/3) implies 

c((i + C)t + 


\\Vp\\l^{l° 


>) - 


< 


c((i+cr ^/^+ (We)’/* .-’/* (‘/* iivvir,j(,.,)’'*). 


Inserting this inequality in the estimate for the vorticity term gives 

Cu 


V 

Ro 


n c puf’ (y/p — 1) dx dr 
! 




+ 


C'z/5/8(We)3/8tV8 ^ ^ 11^2 J |2 \ 

^ II pWlUl^)) 


and, by application of Young inequality, in the end we find 


(58) 


Ro 


0 Jn 


c pu!"^ (y/p — 1) dx dr 




Ro 

+ 


+ 


, 3 r. „^2j 
(Ro) 8/5 ^SWe'l 


2 

LliL^) 


Finally, for the term involving V/^c in (f55l) we can argue in a very similar way. Thanks to the 
bounds of Corollary 13.31 and Lemma l3.5[ where this time we take 5=1/4 and p = 1 (and then 
g = 8/3), we infer 
rt 


V 

Ro 


0 JQ. 


pu-V^c(^- 1 ) 


< 


Cv 

Ro 


II'\/^^IIl7=(l2) IIv^ IIv^IIlRl" 


’) 


< 


Cv 

Ro 


c((l+c)/ + (We)’/* 17711/,,,,)’'' 


Hence, by use of Young inequality as before, it follows the control 
(59) 


Ro 


[ [ PU-'^Uiy/p-l) 
Jo Jn 


< ^C (1 + C) + 


Ro 


^Cv (We) 3/5 ( 8/5 ^^ 31 / 


(Ro) 8/5 


8 We 


I^VII 




Now, we recall equality (f5^ : keeping in mind the definition of (, combining (I58p and (15911 
with (|56p gives us the bound for the general case 1 < 7 < 2; the inequality in the special case 
7 = 2 follows using (f57p instead of (IMp . □ 


Remark A. 10. The rotation term in |13j was dealt with in a slightly different way, exploiting the 
special law of the classical component of the pressure. In the end, one can obtain an analogous 
inequality to the one given here. 
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